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^ ^ ■ The fermionic topological charge of lattice gauge fields, given in terms of 

a spectral flow of the Hermitian Wilson-Dirac operator, or equivalently, as 
the index of Neuberger's lattice Dirac operator, is shown to have analogous 
properties to Liischer's geometrical lattice topological charge. The main new 
result is that it reduces to the continuum topological charge in the classical 
\ continuum limit. (This is sketched here; the full proof will be given in a sequel 

' to this paper.) A potential application of the ideas behind fermionic lattice 

. topological charge to deriving a combinatorial construction of the signature 

\ invariant of a 4-manifold is also discussed. 

o ■ 

O ; PACS. ll.15.Ha, ll.30.Rd, 02.40.-k. - 

'"P' ■ I. Background 

^ . In the continuum, a smooth SU(2) gauge field A{x) = y4^(a:)T°(ix^ on the Euclidean 

^ i hyper-torus is from a mathematical point of view a connection 1-form on a principal 

I SU(2) bundle P over T"^. The bundle is characterised (up to topological equivalence) by 

' an integer Qp and this number can be recovered from the gauge field: 

Qp = Q{A) = -L^^tr(F AF) = ^J d'x e^,,MF,Ax)F^Ux)) (1) 

where F = dA + ^[A,A] is the curvature of A. Thus the topological structure of P is 
encoded in the gauge field A. It is also encoded in the space of zero-modes of the Dirac 
operator = 7'^(9^ + A^): the Atiyah-Singer index theorem [0] gives 

index(^) = Tr(75|^^^^) = QiA) . (2) 

The space of all SU(2) gauge fields on is a disjoint union of components (topological 
sectors) labelled by Q G Z. 
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Now put a lattice (i.e. hyper-cubic cell decomposition) on with lattice spacing 
a. In M. Liischer showed that a lattice gauge field Ufj,{x) , defined on the links [x,x + 
ae^] (e^=unit vector on the positive /^-direction) and taking values in SU(2), also has 
encoded in it the structure of a principal SU(2) bundle P over T^: transition functions 
for P on the overlaps of a collection of regions covering were explicitly constructed 
from U. An integer topological charge for U is then given by Q^'^°{U) = Qp. This is 
often referred to as the geometrical lattice topological charge (to distinguish it from the 
fermionic topological charge discussed below). The construction of P from U is ambiguous 
for certain "exceptional" lattice gauge fields; roughly speaking, these are the fields for 
which it is not possible to canonically write U (p) = exp(r(p)) for each plaquette p , where 
U{p) is the product of the link variables around p. (Example: in the analogous case of 
U(l) gauge fields on the 2-torus 0, U{p) G U(l) can be canonically written as e'^^^^ , 
ir(p) G (— 7r,7r) except for the exceptional fields which have U{p) = —1 for some p.) The 
exceptional fields form a lower- dimensional manifold in the space of lattice gauge fields. 
They can be excluded by restricting to the lattice gauge fields whose plaquette variables 
are sufficiently close to the identity 1 in SU(2); a sufficient condition is 

tr(l-f/(j9)) < 0.03 or \\1-U{p)\\ < 0.015 (3) 

where we are using the matrix norm defined by | |M| p = ^ J2i ioi dxo? matrix M with 
the sum running over the column vectors Mj of M (the normalisation factor ^ is so that 
unitary matrices have norm 1). For the non-exceptional fields the bundle P is uniquely 
determined so Q^^°(U) = Qp is well-defined. Furthermore, the topological structure of P 
is unchanged under continuous deformation of U provided that no exceptional fields are 
encountered . Thus after excising the exceptional fields the space of lattice gauge fields 
acquires a non-trivial topological structure, decomposing into disjoint topological sectors 
which are again labelled by Q G Z. In the continuum Q can take any value in Z , but in 
the lattice setting Q can only take values in a finite subset of Z depending on how fine 
the lattice is f^, §]. When U is the lattice transcript of a smooth continuum field A (and 
U is non-exceptional) the bundle Pu specified by U need not coincide with the bundle Pa 
specified by A in general. However, Liischer showed that Pu and Pa do coincide when 
the lattice is sufficiently fine: he showed that Q^^°{U) —>■ Q{A) in the classical continuum 
limit a ^ [0 . The topology of lattice gauge fields was further elucidated by A. Phillips 
and D. Stone in Q where Q^^°(U) was also generalised to the case of SU(N) gauge 
fields. 

In light of the Index Theorem @) in the continuum, it is natural to ask if there is a 

* For [/ e SU(2) we have tr(l -U) < 2||1 - U\\. 
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lattice Dirac operator D such that 



mdex{Du) = Q''"iU) (4) 

Such a lattice version of the Index Theorem would certainly be mathematically interest- 
ing. It would also be of physical relevance: 't Hooft's solution of the U(l) problem uses 
in a crucial way the connection between the topological structure of A and zero-modes 
of^ implied by the Index Theorem (|^). Also, an alternative "fermionic" description of 
Q^'^°{U) as the index of a lattice Dirac operator could be practically useful for numerical 
work, seeing as the expression for Q^^°{U) is quite complicated and time-consuming to 
implement numerically. Traditional lattice Dirac operators are a bit problematic in this 
context. To avoid doubler zero- modes, operators such as the Wilson-Dirac operator ^ 
include a chiral symmetry-breaking term; this results in the nullspace ker D not being 
invariant under 75 so the zero-modes do not have definite chirality and index D is not de- 
fined. Nevertheless, a lattice version of the index can still be defined from the eigenvectors 
of the Wilson-Dirac operator with low-lying real eigenvalues 0. Essentially the 

same lattice index is obtained Q as minus the spectral flow of the Hermitian Wilson-Dirac 
operator 

HuM = 75(1)^'°" - (5) 

coming from the crossings of the origin by eigenvalues Afc(m) at low-lying values of m.[| 
I (Note that if Xk{rno) = then the corresponding eigenvector ipki^o) is an eigenvector 
for Z}^*'*"" with eigenvalue ^.) This "fermionic" definition of the topological charge of 
a lattice gauge field, which we denote by Qf{U) , can be made more precise by defining 
{U) to be minus the spectral flow of Hu{m) as m varies from to 1; this is justified 
by the fact that for "sufficiently smooth" lattice gauge fields the real eigenvalues of 
jjWiison positive and are localised around — , m = 0, 2, 4, 6, 8 (an analytic explanation 
of this was recently sketched in [^). Note that Q-^{U) has properties analogous to those 
of Q^'^°{U) discussed above: It is defined for all U except those for which Hu{l) has a 
zero-mode; these are exceptional in the sense that they form a lower-dimensional manifold 
in the space of lattice gauge fields. It is clear from the definition in terms of spectral fiow 
that {U) is constant under continuous variations of U provided that no exceptional 
fields are encountered. Thus iU) determines a topological structure in the space of 

'I' The notation in ^ is different; instead of m they have a different parameter K. 

^ In the continuum, taking the 7-matrices to be hermitian so 0a is anti-hcrmitian, crossings of the origin 
by eigenvalues of the hermitian operator HA{m) — 75 — m) only happen at m = and the spectral 
flow from these is — iiidex(j^A) ■ 
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lattice gauge fields in the same way as Q^^°{U) . The exceptional fields for {U) can be 
excluded by imposing a condition on the plaquette variables |]10|, |ll|: 



\l-U{p)\\ < 0.03 



(6) 



Note the similarity with (j^). (Neither (j^) nor are optimal.) 

{U) has subsequently appeared in other guises, as the overlap topological charge 
in |jl2| and more recently as the index of a new lattice Dirac operator: Neuberger's 
Overlap-Dirac operator [|13|, given by 

Hu 



D 



1 



u 



1+75 




(7) 



where Hu = Hjj{l). This operator satisfies |]T4| the Ginsparg-Wilson relation |]T5l Dj^ + 
75D = aD'-f^D , which implies that ker D is invariant under 75 (since Dip = ^ Dl'-fr^ip) = 
{aD'-f^D — 'y^D)ilj = 0) so the zero-modes of D have definite chirality and index D = 

Tr(75|kcrD) 

110,111 



is well-defined, as was first noted in |T^. A formula for the index gives 



mdex{Du) = --Tr(75L>c/) 




(8) 



where the last equality follows from the facts that Tt(^-j^^ is the spectral asymmetry 



of Hu{m) at m = 1 and Hij{m) has symmetric spectrum for m < 

The fermionic lattice topological charge, in its various guises, has been an ongoing 
subject of study since the original works 0, ||, ^ in the mid 1980's; for a selection of 
recent works see, 



e.g. 



21, 19, 22, 23 



24 1 . The emphasis has tended to be on numerical 
investigations though, and there are a number of interesting questions which have been 
partially answered by numerically studies but which currently lack a complete analytical 
resolution. Perhaps the most fundamental of these is the question of whether and Q^*^" 
are equal (at least under suitable conditions on the lattice gauge field, such as (0),(§))-| 
This is equivalent to the question of whether the Lattice Index Theorem (D holds with Du 
being the Overlap-Dirac operator ([^). Another basic question is whether {U) reduces 
to the continuum topological charge Q{A) (= index(j^) ) in the classical continuum limit. 
(This is a necessary condition for equality between and Q^^° since, as discussed above, 
Q3^°{U) — > Q{A) in this limit 0.) This has sometimes been taken for granted in the 



§ This question has received attention in numerical studies, e.g. in 
simpler version of Q3'^°(U) due to P. Woit pst- 



where Qf{U) was compared to a 
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literature; e.g. in it is claimed (without proof) that the eigenvectors of with 
low-lying real eigenvalues will "reduce to" the zero- modes of^ in the classical continuum 
limit. To give a precise mathematical formulation and proof of this statement is not so 
easy though. From a mathematical point of view, showing Qf{U) Q{A) in the classical 
continuum limit is a long-standing open problem, and the main purpose of this article is 
to announce and sketch a proof of this: 

Theorem. For SU(N) gauge fields on the Euclidean 4-torus, iU) Q{A) in the classical 
continuum limit. Equivalently, index(D[/) — > index(^A) in this limit, where Du is the 
Overlap-Dirac operator (|^) above. 

This is the fermionic analogue of Liischer's result Q^^°{U) —>■ Q{A). The deeper 
question of whether and Q^'^° are equal (at least when U{p) satisfies a bound of the 
form (^,(§1)) is a challenging mathematical problem which we will not attempt here. 

II. Classical continuum limit of the fermionic lattice topological charge 



The proof of the theorem above, which we sketch in the following, grew out of 



suggestions by Martin Liischer . The full details will be given in a forthcoming paper 
P7[ . (An alternative argument for a restricted class of topologically non-trivial fields in 
a slightly different setting]]] was previously given in [^, v4].) The general idea is to start 
with the last equality in (H), 



Q^(^) = -^Tr(^), (9) 
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carry out a certain power series expansion of the inverse square root, and then evaluate 
the trace in an explicit basis for the spinor fields. But we do this in a slightly indirect 
way: First, we use the locality result of |^ for the Overlap-Dirac operator to derive a 
relation Q^{U) = Qf {U)^'^p+^^ + O^e-^'") (with c> 0) where Q-^(f/)(2p+i) is the fermionic 
topological charge in a setting in which an infinite volume limit p oo can be taken. 
Then lima^o {U) = lima^o hnip^oo Q-^(f/)*^^^+^^ and the latter quantity is easier to 
evaluate because the sum resulting from the trace in (|^) becomes a tractable integral in 
the p — oo limit. 



II-l. Preliminaries 

The 4-torus T'^ is taken to be [—^L, \L\/ ^ x ■ ■ ■ x [—\L, \L\/ (where ~ means 
identify endpoints). Then a continuum SU(2) gauge field v4^(x) on can be viewed as 

^ which can be reduced to the present setting Q 
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a gauge field on R'' satisfying 

Af,{x + Le^) = z/)A^(x)^](x, z/)"^ + u)'^ (10) 

wfiere ^(x, z/) , u = 1,2,3,4, are tlie SU(2)-vaIued monodromy fields which specify the 
principal SU(2) bundle over T^.^ Now put a hyper-cubic lattice on with 2N sites along 
each edge, so the lattice spacing is a = L/2N. This extends to a hyper-cubic lattice on 
with sites aZ^. The lattice transcript of A , 

t/^j(x) = Texp(^ J aAfj_{x + tae^) dt^ (11) 

(T=t-ordering and for simplicity the coupling constant has been set to unity) satisfies 

U^{x + Ley) = fi(x, u)Uf,{x)n{x + ae^, z/)"^ . (12) 

The finite-dimensional complex vector-space C of lattice spinor fields on consists of the 
spinor fields ip{x) , x G aZ*^ , satisfying 

ip{x + Le^) = Q{x,h')ip{x) . (13) 

The covariant finite difference operators ^V^ , given by 

V+^ix) = Uf,{x)^{x + acf,) - tPix) (14) 
V~ip{x) = i){x) -U^{x - ae^Y^ipiyX - ae^) , (15) 



preserve (0) and are therefore well-defined on C , as is the Wilson-Dirac operator, given 

byQ 

j^Wilson ^ ^ |a(^Af;) (16) 

where ^ = ^I]^7^V^ = ^Z]/i7^|(V^ + V~) is the naive lattice Dirac operator 
(the 7'^'s are taken to be hermitian so is anti-hermitian), = ^I]^V~ — = 
^S/i(V^)*V^ = ^Z]^(V~)*V~ is the lattice Laplace operator (hermitian, positive). 
Likewise Hjj = 75(1^^*'*°" — i) is well-defined on C , and so is the Overlap-Dirac operator 

II These also satisfy a cocycle condition which ensures that Af^{x + Le^, + Lcp) is unambiguous. It is 
always possible to make a gauge transformation so that fl(x, ly) ~ I for = 1, 2, 3 and Q{x, 4) is periodic 
in xi,X2,X3. Then for fixed X4 il{x, 4) determines a map — + SU(2). The degree of this map (which 
is independent of X4 since il(x,4) depends smoothly on 2:4) is precisely the integer Qp specifying the 
SU(2) bundle P over T"* , and is therefore also the topological charge Q{A) = Qp o{ A. 

** For simplicity we have taken the Wilson parameter to be r = 1. 
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Djj = ^(1 + '~i5Hu\Hjj\^^) provided that Hu has no zero-modes. Hu is guaranteed not to 
have zero-modes when a is sufficiently small. Indeed, 1 — U{px,iJ.,i) = —c^^F^y^x) + 0{a?) 
vanishes uniformly for a — > , implying that (H) is satisfied for sufficiently small a , which 



in turn implies that [TT| 

> b > (17) 

(This holds, e.g., with 6 = 0.1 but we will not need the explicit value in the following.) 
We henceforth assume a to be small enough that (^ - and thereby (|T^) - holds. 



II-2. Passing to a setting with an infinite volume limit 

A general linear operator O on lattice spinor fields corresponds to a kernel function 
0{x,y) via 0%Ij{x) = a'^Y^yO{x,y)ip{x). By (|) we have 

Qf(U)=a^Y.Mx) (18) 

where 

qu{x) = -^tii'y^Duix^x)) (19) 

and the summation is over 

T = {x = an\nf, = -N, -A^ + 1, . . . , - 1} . (20) 
Now, for arbitrary whole number p , set 

= Q{x + pLe„, z/)fi(x + (p — l)Le,y, u) ■ ■ ■ Q{x + Le,y, i')Q{x, v) (21) 
and define Cp to be the space of lattice spinor fields '?/'(x) , x G aTJ^ , satisfying 

%l){x + (p + l)Le^) = fi(P)(x, iy)tlj{x) . (22) 



Note that (p!3| ) implies (|2^), so C is contained in Cp for all p. Note also that the covariant 
finite difference operators (|I^ ) -([T5| ) preserve (|2^) and are therefore well-defined operators 



on Cp ; it follows that the Overlap-Dirac operator is well-defined as an operator on Cp. We 

e fo 



denote the Overlap-Dirac operator on Cp by D^j^ in the following, with Du denoting the 



operator on C. The fact that Du is the restriction of D j to C C Cp for all p implies 

Du^{x) = DS'"''^^{x) = a' Y: DS'''-'\x,y)^{y) for ^ G C (23) 
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where 

Tsp+i = {x = an\n^ = -{2p+l)N,-{2p+l)N +l,...,{2p + l)N -1} 

= {x + Lm I X G r , = —p, —p + 1, . . . ,p} . (24) 
From this it is straightforward to derive using (|T^, that for p > 1 the norm of 

R^^'^'\x,y) := Du{x,y) - Df^'\x,y) (25) 
has a bound 

\\B!i''^'\x,y)\\ < E \\Df'-'\x,y + Lm)\\. (26) 

|m^|G{l,2,...,p} 

The locahty resuh in now gives 

||aDg^+^)(x,a; + Lm)|| < ^ exp(-C2- ^ |m^|) (27) 

where Ci and C2 > are constants independent of a, p and U. It follows that 
||ai?g^+^)(x,x)|| < 2^ E ^exp(-C2-E^. 

m^G{l,2,...,p} M 

a* „ Ji/2 ^ a 



'1/2 



Set 



qS'^'\x) = -%{i,D^J'^'\x,x)) (29) 
then substituting (^5]) in ([T9|) and taking account of (|28|) we see that the norm of 

r(,^^+^)(x):=gf;(x)-g(,^^+^)(x) (30) 
has a bound of the form 

\r'u'"^^\x)\ < Oie-'/") OO (31) 
where 0{e~'^^"') is independent of U and p. It follows that 

lim qu(x) = lim lim qlf^~^^\x) (x G T) (32) 

lim g-^(f/) = lim lim a^yg^'*'+')(x) (33) 

xGr 

The infinite volume limit p — > oo in these expressions will facilitate their evaluation, as 
we will see in the following subsection. 
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II-3. Evaluation of the classical continuum limit 

We now exploit the locality of the Overlap-Dirac operator in the gauge field |]TU to 
replace or rather, its lattice transcript f/^(x), in (p2D -(^) by the lattice transcript 

Ufi{x) of a gauge field 



X 



which coincides with Af^(x] 



in a neighbourhood of [—^L, |L]^ 
but which vanishes outside a bounded region in R^. Specifically, choose a smooth function 
X{x) on which equals 1 in a neighbourhood of [—^L, \LY and vanishes outside a region 



contained in [—^L, |L] 



and set A^^^x) 



Kx)A^ 



\x] 



For each p > 1 we take U to be the 



lattice transcript of A in [— ^^^L, ^^^L]^ and extend U to the rest of the lattice on R^ by 
requiring that U^{x) be periodic in all directions with period (2p+ 1)L. Then the Overlap- 
Dirac operator with lattice gauge field f7 is a well-defined operator D-^^^'^ on the space 

C2P+1 of lattice spinor fields satisfying the periodicity condition ip{x + {2p+ l)Le,y) = ip^x) , 
u = 1,2, 3, 4. A version of the locality result of [|10] leads to p7[] 



\aDS^+'\x,x) 



aD^y^^Hx,x) 



< 



where c > is a constant and 0(e ^Z'*) is independent of x 
l2D-f53D are unchanged if qij^^^\x) is replaced by 
1 



xeV p>l (34) 
U . It follows that 



U . 



(2p+l)/ X 
9~ W 



-tr(754'^+'^(x,x)) 



There are several reasons for making the replacement A A , U - 
justification in of the steps that follow uses the fact that A^[ 



X 



(35) 

U . (i) The rigorous 
) vanishes outside a 



bounded region. (This is not the case in general for A^[x) which can diverge in the 
|x| — *• cxo limit.) (ii) It leads to Cp being replaced by a space Cp of periodic spinor fields, 
thereby allowing the trace in (|36|) below to be evaluated in a "plane wave" basis. 

Substituting the expression for the Overlap-Dirac operator in ( pSj) we find (cf (|D) 



(2p+l)/ N 



1 



-tr 



u 



{x,x] 



where 



1 

— tr 
2 



-1 
2^ 



}(x,: 



(36) 



X; 



u 



( T~,Wilson 1\ 



(37) 
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and the operator is defined by {Sx4'){y) = ip{x)5xy The strategy now is to carry out a 
power series expansion of the inverse square root in (0). A calculation gives 

^h^u = Lu + V^ (38) 

where 

L = _V' + i(A-2)2 (39) 

V = -i[7^71[V,,V.]-i^,A]. (40) 

is a linear combination of commutators of the V^'s. As pointed out in [|1^], the norms 
of these commutators are bounded by maxp||l — t/(p)||. In the present case, this together 
with the expansion 

f/(p.,^,.) = 1 + a^F^uix) + Oia') (41) 
show that 

IIV^II ~ Oia') (42) 



As in { \L7\j we have a bound X~X^ > b > when a is sufficiently small; furthermore, 
due to (^) we can assume that \\V\\ < b/2. Then from ( p8[ ) we get a lower bound 
Lfj > 6/2 > on the positive hermitian operator L^. It follows that L is invertible and 
||L~^\^|| < 1 when a is sufficiently small. The inverse square root in (PBf ) can then be 
expanded as 



oo 



da 1 



'•'^ da / 1 



l"^ da / i w 
i-oo V V L + a2^ V I 



+ (L + a2)-iV 



J 1 

5:(-l)'=((L + a^)-V)^ (43) 

fc=0 



-oo TT L + O"^ 



Note that the 7-matrices in ( ^3]) are all contained in V. Since the trace of 75 times a 
product of less than 4 7-matrices vanishes, the terms with A; = and = 1 in (|4^) give 
vanishing contributions to (PB]). On the other hand, by (^2]) the terms with k > 3 are 
O(y^) ~ O(a^) , hence the contribution of these in ( p6D is ~ O(a^). (This is rigorously 
established in using the presence of 6x in (^) together with the fact that A^{x) has 
compact support.) Thus the only relevant contribution to (|36| ) from the expansion ( ^3] ) 
comes from the k = 2 term: 

°o drr 1 1 

- 2^2 ^^^^ (44) 

vr (L + 0-2)2 + ^ ^ 
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By making a power series expansion of + cr^) and using [L, V] ~ 0{a) we find in [27 
that, modulo an 0{a) term, the contribution of (|^) in (|36D is the same as that of 



da 



3 vr (L + a2)3 
We hereby see that (BH) reduces to 



J — ( 



da 



TT fl + a2)3 



y 2^-5/2 



-3 
160^ 



Tr(75X^V|L~^/'5.) + 0(a) 



(45) 



(46) 



Now, from (|TT]) with A ^ A we see that the a-dependence of U enters through a product 
aA. It follows that 



L^ = Li + 0{a) 



(47) 



Since V ~ O^a"^) we can replace Xi and L^^ — Li in ( ^6[ ) at the expense of another 

0{a) term. With the expressions (|37| ) and ( ^9|) for X and L we find 



q^~'^'\x) = T7r^Tr(754Vg(yi + |(Ai - 2))(-V? + |(Ai - 2))-^/^) + 0(a) (48) 



The rigorous derivation of this in |^ again uses the fact that v4^(x) has compact support. 
The trace in (?8) can now be evaluated using the "plane wave" orthonormal basis {(pk} 
for the lattice scalar fields with periodicity 0(a; + (2p + l)Leu) = 0(x) (i.e. the fields in 
the scalar version of C2p+i), given by 



Ak-x 



M ={{2p+l)2NY 



(49) 



a(2p + 

Note that the volume per k is 

27r x4 (27r)^ 



-|^{ - {2p + l)iV, -(2p + l)iV + 1, . . . , (2p + l)iV - 1} (50) 



A^A;= ( 



a(2p + 1)2X 



aW 



(51) 



Using the plane wave basis, the trace in ( ^8]) can now be evaluated as in [^, v4], leading 
to 



9: 



(2p+l) I 
!7 



X 



) = q^{x) J2 a^^^k i{ak) + 0(a) 



(52) 
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where the summation region for k is (0), A'^k is given by (^Tj), 

. . 1 



and 



32 



3 n,.cosir^(-l + - cos/c^) - E 



^^ cos fc„ 



Eu sin^ + (-1 + ^ cos A;^))2 



5/2 



(53) 



(54) 



and we have exploited the fact that A^(a;) = Af^{x) for a; G F. Changing summation 
variable from k to k = ak in (K^) and taking the p oo limit gives 



lim gl'^+')(x) = qA{x) I d^kl{k) + 0(a) 



(55) 



The integral over k in this expression was evaluated in pSf, and independently in p9| , 
and was found to be 1. Recalling (p^)-(|33|) we finally get 



lim qu{x) = lim lim q\j 



(2p+l), 



a->0 



a-*0 p-*oo 



x) = lim lim q~^^^\x) 



[x e r) (56) 



and 



lim Q^(f/) = lim lim a^yg}f+^^(x) 



lima^y (qA(x) + 0(a)) = [ d^x qA(x) = Q(A) 



(57) 



This completes the sketch of the proof of the theorem. The latter part of the derivation 
has similarities with, and was inspired by, the calculation of the classical continuum limit 



of the axial anomaly for Wilson-Dirac fermions done many years ago by W. Kerler [30 



and E. Seller and I. O. Stamatescu ||3T[| . 

The Overlap-Dirac operator determines a lattice Dirac fermion action ipDuip |]T3 



which has an exact symmetry under a new kind of lattice chiral transformations [|T^ . The 
corresponding axial anomaly (=the infinitesimal chiral jacobian) was found in [|l^ to be 
(in the flavour singlet case) 

Au(,x) = 2iqu{x) 

Its classical continuum limit was calculated in a perturbative setting in and subse- 
quently in a non-perturbative setting in vl],||T8|,|29[. The perturbative setting was 
further studied in [^. However, these calculations are all problematic in the case where 
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the continuum gauge field A^{x) is topologically non-trivial (e.g. A^{x) does not have a 
well-defined Fourier expansion in this case.) This case is covered by the arguments above 
though: we have found that 

quix) = qA{x) +0{a) 

in complete generality. As we have seen, the locality of the Overlap-Dirac operator (which 
was not used in the previous references) is a key ingredient in the derivation of this result. 

The arguments and results above, which are for SU(N) gauge fields on the 4-torus, 
can be easily extended to U(l) fields on the 2-torus, thereby providing analytic verification 



of results from numerical studies in [0 where the fermionic topological charge and axial 
anomaly were seen to reduce to the correct continuum quantities for topologically non- 
trivial gauge fields in this setting. 

III. Index of lattice Dirac operators and combinatorial approach to topological 
invariants 

In this section a potential mathematical application of the ideas behind fermionic 
lattice topological charge to the construction of certain topological invariants of manifolds 
is discussed. Consider Kahler-Dirac spinor fields ^ with the "spacetime" being 
an arbitrary smooth oriented riemannian 4-dimensional manifold M. In the continuum 
the fields are the differential forms on M, the space of which we denote by Q{M) = 
®^^(,nP{M). (Locally, a p-form 00 G nP{M) is of the form ujf,,,„i,^{x)dx^'^ A ■ ■ ■ A dx''".) 
The Dirac operator on Q{M) is 

D = d + d* (58) 

where dp : ^^(M) QP'^^IM) is the exterior derivative. The standard way to define 
the chirality operator Fs (analogue of 75) is F5 = (—1)^ on Qp{M). This has the chi- 
rality properties (Fs)^ = 1 and T^D = — -DF5 , and it is not difficult to show that the 



corresponding index of D is the Euler characteristic xi^) of M (see, e.g., |]36| , Ch.4]) 



A discretisation of Kahler-Dirac theory can be obtained via a triangulation i^' of M 



35| . In the discrete theory the fields are the cochains of K (=the functions on the 
oriented simplexes of K), the space of which we denote by C{K) = (Bp^oC'^iK). The 
analogue of d is the coboundary operator d^ : Cp{K) Cp~^^{K) , the Dirac operator 
is = d^ + {d^y and F5 = (—1)^ on Cp{K). Using the de Rham theorem it can be 
shown that index = index D = x(M) so index is a combinatorial construction of 
the topological invariant x{M). This is nothing new though; in fact we have just recovered 
Euler 's original combinatorial construction of x{^)- Now, there is another way to define 
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Fs in the continuum, namely by = (—1)'^* where * : Qp{M) — > f2^~''(M) is the Hodge 
star operator and s = (— on f2P(M) ; this also satisfies the chirality properties 
(Fs)^ = 1 and T^D = — -DF5 and the corresponding index is again a topological invariant 



of M, namely the Hirzebruch signature a{M) (see, e.g., |]36| , Ch.4]). Unlike xi^) there 
is no known combinatorial construction of a{M) - this is an interesting open problem in 
mathematics. The problem is refiected in the fact that there is no natural discretisation 
of * , and thereby of F5 , in this setting. However, a natural discretisation is obtained 
after introducing K : the cell decomposition of M dual to K. We can then consider the 
doubled discretisation f2(M) C{K) © C{K) with 



where *^ : Cp{K) C^~p(K) and *f : Ci(K) C^-'i{K) are the duality operators 
(see, e.g., where topological quantities were exactly reproduced in a similar discrete 
setting). The index of can be seen to vanish though; this is essentially due to the 
fact that rf^ is skew-diagonal in (|59|) . This is reminiscent of the vanishing of the index of 
the usual naive lattice Dirac operator. The necessity of introducing C{K) is reminiscent 
of the appearance of fermion doubling in the usual lattice theory. Therefore, in light of 
the preceding sections, it may be possible to obtain a combinatorial construction of a{M) 
from the spectral fiow of a "Hermitian Wilson-Dirac operator" = {Dwuson ~ 

m/(a)).[^ The main problem is to find a suitable "Wilson term" W^^ on C{K) © C{K) 
for the 'Wilson-Dirac operator" D^^fi^^^^ = iD^^ + aW^^ (where a is now the mesh size 
of K), and a suitable function /(a) in . (In the usual lattice setting /(a) = 1/a.) 
One idea is to embed C{K) ®C{K) into C{BK) where i?i^=the barycentric subdivision 

of and then try to construct W^^ from the discrete Laplace operator A^^^ on C {BK) . 

Finally, we remark that after "twisting" by a fiat gauge field A, the signature (j^(M) 
of a 4-manifold M with boundary is closely related to the Atiyah-Patodi-Singer rho 



invariant p(A^, a) of a 3-manifold together with a representation a : 7ri(A^) 0{n) |38 
A combinatorial construction of the signature invariant (Ja{M) would in all likelihood lead 
to a combinatorial construction of the rho invariant as well. 
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